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Abstract 

Recently there has been some progress in modehng the statistical properties of turbulent flows 
using simple superstatistical models. Here we briefly review the concept of superstatistics in 
turbulence. In particular, we discuss a superstatistical extension of the Sawford model and 
compare with experimental data. 



Turbulence is a spatio-temporal chaotic dynamics generated by the Navier-Stokes equation 

il= -{vV)v + uAv + F. (1) 

In the past 5 years there has been some experimental progress in Lagrangian turbulence mea- 
surements, i.e. tracking single tracer particles in the turbulent flow. Due to the measurements 
of the Bodenschatz [1, 2, 3] and Pinton groups [4, 5] we now have a better view of what the 
statistics of a single test particle in a turbulent flow looks like. The recent measurements have 
shown that the acceleration a as well as velocity difference u = v{t + t) — v{t) on short time 
scales r exhibits strongly non-Gaussian behavior. This is true for both, single components as 
well as the absolute value of a and H. Moreover, there are correlations between the various 
components of a, as well as between velocity and acceleration. The corresponding joint prob- 
abilities do not factorize. Finally, the correlation functions of the absolute value |o| and |u| 
decay rather slowly. 

How can we understand all this by simple stochastic models? There is a recent class of 
models that are pretty successful in explaining all these statistical properties of Lagrangian 
turbulence (as well as of other turbulent systems, such as Eulerian turbulence [6, 7, 8], atmo- 
spheric turbulence [9, 10, 11] and defect turbulence [12]). These are turbulence models based 
on superstatistics [13]. Superstatistics is a concept from nonequilibrium statistical mechanics, 
in short it means a 'statistics of statistics', one given by ordinary Boltzmann factors and an- 
other one given by fluctuations of an intensive parameter, e.g. the inverse temperature, or the 
energy dissipation, or a local variance. While the idea of fluctuating intensive parameters is 
certainly not new, it is the application to spatio-temporally chaotic systems such as turbulent 
flow that makes the concept interesting. The first turbulence model of this kind was intro- 
duced in [14], in the meantime the idea has been further refined and extended [15, 16, 3, 8]. 
The basic idea is to generate a superposition of two statistics, in short a 'superstatistics', 
by stochastic differential equations whose parameters fluctuate on a relatively large spatio- 
temporal scale. In Lagrangian turbulence, this large time scale can be understood by the fact 
that the particle is trapped in vortex tubes for quite a while [3]. Superstatistical turbulence 
models reproduce all the experimental data quite well. An example is shown in Fig. 1. The 
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Figure 1: Probability density of an acceleration component of a tracer particle as measured 
by Bodcnschatz ct al. [1, 2]. The solid line is a theoretical prediction based on lognormal 
superstatistics (s^ = 3) [16]. 



theoretical prediction which fits the data perfectly is given by 

^2 



with n = 62*^ and only one fitting parameter, = 3.0. A similar formula as eq. (2) was 
already considered in [17], though without a dynamical interpretation in terms of a stochastic 
differential equation with fluctuating parameters. 

The key ingredient of superstatistical models is to start from a known model generating 
Gaussian behaviour, and extend it to a superstatistical version exhibiting 'fat tails'. In 
general, in these types of models one has for some dynamical variable a the stationary long- 
term density 

p{a) = ^me-h^'^'dp, (3) 

where is some suitable probability density of a fluctuating parameter /3. The function 
f{(3) fixes the type of superstatistics under consideration. In particular, it is responsible for 
the shape of the tails [18]. Note the mixing of two statistics, that of a and that of p. 

In Lagrangian turbulence, one may first start from a Gaussian turbulence model, the 
Sawford model [19, 20]. This model considers the joint stochastic process {a{t),v{t), x(t)) of 
an arbitrary component of acceleration, velocity and position of a Lagrangian test particle, 
and assumes that they obey the stochastic differential equation 



l2al{Tj;'+t^')TzH^' L{t) (4) 
V = a (5) 
X = V, (6) 
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L{t): Gaussian white noise 

Tl and t^f. two time scales, with Tl » tfj, 

Tl = lolliC^t) 

e: average energy dissipation 

Co,ao: Lagrangian structure function constants 

cr^ variance of the velocity distribution 

Rx = a/T5<7^ / Taylor scale Reynolds number. 
For our purposes it is sufficient to consider the limit Tl ^ oo, which is a good approximation 
for large Reynolds numbers. In that limit the Sawford model reduces to just a linear Langevin 
equation 

a = -7a + crL(i) (7) 

with 

7 = (8) 
^3/2 

a = ^ly-'/h. (9) 
2ao ^ ^ 

Note that this is a Langevin equation for the acceleration, not for the velocity, in marked 
contrast to ordinary Brownian motion. 

Unfortunately, the Sawford model predicts Gaussian stationary distributions for a, and is 
thus at variance with the recent measurements. So how can we save this model? 

As said before, the idea is to generalize the Sawford model with constant parameters to 
a superstatistical version. To construct a superstatistical extension of Sawford model, we 
replace in the above equations the constant energy dissipation e by a fluctuating one. One 
formally defines a variance parameter [16] 

27 4ao 1/2 J_ . . 

^' C2 £3/2' ^^^^ 

where e fluctuates. Now, if (5 varies on a large spatio-temporal scale, and is distributed with 
the distribution f{0), one ends up with eq. (3) describing the long-term marginal distribution 
of the superstatistical dynamics (7). This is basically the type of model introduced in [14], 
there with chosen to be a ^^-distribution. Models based on x^.g^pgi-g^atig^icg yield good 
results for atmospheric turbulence [9, 10], and ultimately lead to Tsallis statistics [21]. On 
the other hand, for laboratory turbulence experiments one usually obtains better agreement 
with experimental data if /(/?) is a lognormal distribution. In view of eq. (10) this is clearly 
motivated by Kolmogorov's ideas of a lognormally distributed e [22]. 

Superstatistical models are not restricted to Lagrangian turbulence but can be also for- 
mulated for Eulerian turbulence [7, 14]. Fig. 2 shows that also here one obtains excellent 
agreement with experimental data: Probability densities p{u) of longitudinal velocity differ- 
ences u are well fitted by lognormal superstatistics on all scales. The parameter varies 
with the scale. In fact, not only the distribution p{u) but also the distribution /(/?) can 
be directly measured in experiments [8], and the two can be consistently connected via the 
superstatistics formalism. Jung and Swinney [8] have also experimentally confirmed a simple 
scaling relation between /? and the fluctuating energy dissipation e. 
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Figure 2: Experimentally measured histogram of velocity differences in a Taylor-Couette 
experiment [6], and comparison with a superstatistical prediction (s^ = 0.28). 

It should be noted that if we know the probability densities p{u) analytically, as well as 
the dependence of the parameter on the scale r, we can also calculate moments of velocity 
differences and thus determine scaling exponents (m defined by 

Many different models of can be constructed in such a way [16, 23, 24]. For further 
stochastic models, see e.g. [25, 26, 27]. 
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